On the current gap of single-electron transistors 
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Effects of the single-electron tunneling and the Coulomb blockade in a cluster structure (the 
molecular transistor) are investigated theoretically. In the framework of the particle-in-a-box model 
for the spherical and disk-shaped gold clusters, the electron spectrum, the temperature dependence 
of the chemical potential and the residual charge are calculated. We show that the residual charge 
is equal to the non-integer value of elementary charge e and depends on the cluster's shape. The 
equations for the analysis of the current-voltage characteristic are used under conservation condition 
for the total energy of the structure taking into account the contact potential difference. Restrictions 
associated with the Coulomb instability of a cluster are introduced into the theory in a the simple 
way. It is shown that the critical charge of the cluster in the open electron system is close to the 
residual charge. For single-electron transistors based on small gold clusters the current gap and its 
voltage asymmetry are computed. We demonstrate that the current gap exhibits non-monotonic 
size dependences, which are related to the quantization of the electron spectrum and the Coulomb 
blockade. 

PACS numbers: 72.20.-i, 73.22.-f, 73.23.-b 



-a 
c 

o 
o 



> 
m 
m 

(N 

in 
o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

Metal granules, which are weakly coupled via tunnel 
barriers to electron reservoirs, are of considerable interest 
in physics of low-dimensional systems (see reviews 
and references therein). 

The tunneling current flowing through two massive 
electrodes can be controlled, if a cluster is placed between 
them. At first sight, the probability of electron tunnel- 
ing (and consequently a value of the current) should be 
much greater in the presence of a granule between the 
reservoirs, than in the case of its absence. However, an 
opposite behavior was observer in experiments for the 
spherical-like [IE] IS El an d disk-shaped 0,13 clusters. 
Measured I — V characteristics have a plateau of the zero 
current (a current gap). 

In Refs. 0, d, @ the structure with two tunnel junc- 
tions (Fig. ^) was represented by a thick Au (111) film 
covered by the d e ~ 10 A thick dielectric layer (with the 
dielectric constant e « 3), on which the small spherical- 
like gold clusters were organized. The tungsten tip (with 
small curvature of surface) of STM microscope was cov- 
ered by a gold film with 10 3 A thickness. Therefore 
we can consider all the three electrodes (two of them 
with a flat surface) as being gold. Using the circuit ap- 
proach of Ref. 9], the capacitances, tunnel resistances 
and a "residual" (fractional) charge Qq of granules were 
extracted as fitting parameters from the measured de- 
pendence /(V"). 

Earlier, a similar STM (Pt/Ir -tip) measurement for 
such a structure (d e ~ 1.4 nm, e ~ 2.7) based on a gold is- 
lands with monatomic height H « 0.25 nm (disk-shaped 
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FIG. 1: The energy diagram for the structure Au/Au4o/Au 
before mpJication_of_ voltage; ei = e, £2 = £3 = 1 in experi- 
ments I 



clusters) was carried out in the works 0] • 

The experiments demonstrated the following features 
of the I(V) behavior: 

1. The gap width of the zero conductance is approx- 
imately proportional to the inverse radius of the disks 
(Fig. 4 in Ref. 0) and spheres (Figs. 1(c) and 2(a) in 
Ref. 0). This does not allow one to establish unequivo- 
cally classical or quantum origin of the gap. On the other 
hand, out of the current gap, the steps of the staircase 
are clearly visible (Fig. 3 in Ref. [jj and Fig. 1(b) in 
Ref. 0). 

2. For a disk, the gap width varies non-monotonically 
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with alteration of the collector-cluster distance under the 
fixed emitter-cluster one (Fig. 3 in Ref. (g)- 

The possibility of a fractional charge at tunneling 
structures was discussed in the Refs. [l], ^j. In the 
percolation systems it is supposed that the charge Qo at 
the each granule has a soliton origin. The value of this 
charge was calculated numerically in Refs. [TllIT^. 

However, not much attention was paid to this problem, 
as well as the current gap origin on the I — V curves and 
its asymmetry. Perhaps, this problem is related to the 
fractional quantization (or fractional statistics) , when the 
decoupling of the spin and the electron quantum numbers 
of a charge is important. 

The aim of this work is the computation of the current- 
voltage characteristic of the three-electrode structure, 
whose central electrode is a metal cluster with different 
sizes and shapes. 

The outline of this paper is as follows. In Sec. II we for- 
mulate the problem of the influence of the well-separated 
energy levels and charging energy on the resonant tunnel- 
ing through a metal cluster, which is weakly coupled to 
two massive electrodes. Our main assumption is that the 
thermal energy exceeds the width of transmission reso- 
nance. The present approach describes the case of strong 
inelastic scattering of electrons in the clusters that cor- 
responds to the full thermalization regime. We calcu- 
late electron spectrum for the wells of various shape, 
chemical electron potential and residual charge. This 
charge is a result of equilibrium with electron sea. We 
use the density-functional approach for determination of 
the total energy of metal cluster in an external electric 
field. Next, we estimate the critical surplus charge of 
the spheroidal cluster, which leads to the Coulomb in- 
stability for a nonzero applied voltage. In Sec. Ill, we 
consider electron transport through a metallic quantum 
dot that can be described by a means of master equation. 
The contact potential difference is taken into account. In 
Sec. IV we perform the numerical analysis of the I — V 
curves for various set of parameters. The current jumps 
are investigated in details for the tunneling structures 
based on the magic and non-magic clusters. The model 
allows one to determine the size of current gap and its 
asymmetry on a voltage. 



II. PRELIMINARY ANALYSIS AND 
FORMULATION OF THE PROBLEM 

We consider spherical gold clusters whose radii vary in 
the range R ~ 7 - 14 A, R = N 1/3 r s => N ~ 100 - 
600 (r s — 3.01 ciq is the electron density parameter, ao 
is the Bohr radius). Similarly for disks of monatomic 
thickness: R ~ 5-42.5 A^ N - 14- 10 3 . We introduce 
the characteristic charging energy Eq = e 2 /C, where C 
is a cluster capacitance 13] . For spheres and disks we 
obtain E G ~ 1.82 - 1.06 and 3.2 - 0.42 eV, respectively. 
Temperatures of structures are T < 30 K. 

Let's determine the electron spectrum in spherical and 



cylindrical wells (see Appendix A). For the above men- 
tioned sizes, calculation in both cases gives the close val- 
ues of the spectrum discreteness, Ae p « 1.2 — 0.3 eV, 
nearby to highest occupied level e HO at T = 0. Thus, for 
the whole range of R in experiments 0, ID, El we have 
to deal with a set of open 0D systems (quantum dots) . 
The resulting inequality, 

E c « Ae p > k B T, (1) 

corresponds, apparently, to two coexisting structures at 
I — V curves: effects of the spectrum quantization and 
the Coulomb blockade. The current-voltage character- 
istic should represent a superposition quantum staircase 
with a step ~ Ae p /e and the classical Coulomb one of the 
electrostatic nature with a step ~ En le a long a voltage. 
However, detailed measurements 011 till Eifperformed 
to date do not yield an unequivocal conclusion about the 
effect of electron quantization levels upon the /(V). 

In our opinion, the discreteness of the spectrum actu- 
ally determines the zero conductance gap of the I — V 
curves observed in 0, IE El 0> 13 • Let's consider the prob- 
lem step by step. 

A. Structure in the absence of voltage 

The left and right electrodes (emitter and collector) 
represent the electron reservoirs. Each reservoir is taken 
to be in thermal equilibrium. A continuum of states is 
assumed in reservoirs, occupied according to the Fermi- 
Dirac distribution 

f(e e > c + W^ c ) = {1 + exp[( £ - c + W^n/kBT}}- 1 , (2) 

where Wo > is the work function for a semi-infinite 
metal. In all cases energy Z7 < £ < is counted off from 
the vacuum level, Uq < is the position of conductivity 
band of a semi-infinite metal. 

The electron chemical potential fi of a cluster in a 
quantum case can be defined by the normalization con- 
dition 

oo 

£/(e p -A0=tf, (3) 
p=i 

where sum runs over all one-electron states, N is the 
total number of thermalized conduction electrons in a 
cluster (taking into consideration the surplus and lacking 
electrons), and 

f(e p - /i) = {1 + exp[(e p - Ai )/fc B T]}- 1 . (4) 

If the electron spectrum is known, from Eq. @ it is 
possible to calculate fi of Aujv , where No is the num- 
ber of the conduction electrons of neutral clusters. The 
Fermi level of non-magic clusters coincides with a real 
level in a cluster. For the magic ones it lies between the 
energy terms. Figure |21 depicts the chemical potential 
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FIG. 2: Temperature dependence of chemical potentials of 
the spherical gold non-magic AU39 and magic AU40 and Aui96 
clusters. 

of some spherical clusters as a function of temperature. 
Predictably, the dependence is slack and is completely 
determined by the level hierarchy in dots, and also by 
the number of electrons. Calculations show, that the 
temperature gradient of chemical potential can be both 
the positive and the negative, and at some temperatures 
it can change a sign. Similar behavior p,(No,T) for magic 
clusters Na^v have been reported in Ref. |15j . 

A contact potential difference appears between a clus- 
ter and electrodes (Fig. 

5<t>=(W +(i)/e. (5) 

An equilibrium is reached by the noticeable charging of a 
cluster since it capacitance is finite. If \fi\ < Wq, a cluster 
is charged positively by a charge Q = —e(N' — N ) > 0, 
where N' is determined by the solution of the Eq. © 
with replacement /i — > — Wq for a spectrum {e p } shifted 
on —e54>, according to the Koopmans' theorem |13(. In 
similar theories, the presence of a contact potential dif- 
ference were neglected. 
Thus, we have 

Qo = C5<t>. (6) 

This expression differs from the definition of Qo in the 
orthodox theory (Eq. (29 b) in Ref. [T^l. 

A quasi-classical approximation, Wo + n(R) — fii/R, 
\i\ = 1.9 eVxao pj, gives Qo — +0.07 e. In a quantum 
case, it is necessary to make a replacement fi — > e H °. 
Otherwise, when > Wo the cluster is charged nega- 
tively [p — > e LU ). 

The non-integer value of Qo ( m terms of the elemen- 
tary charge) is due to the fact that the electron wave 
functions are not well localized. Therefore, electrons can- 
not be interpreted as classical particles, and a fraction of 
an electron (and its charge) can be found in the other 
electrode 0,113. 

Value Qo ~ +0.5 e better than the others corresponds 
to the Kuzmin's and Likharev's experiment (Fig. 2(b) 



in Ref. 0) m which the I — V was measured for the 
structure of two electrodes (alloy Pb/Au) and granule In 
of radius R « 100 nm separated by the oxide films. It is 
interesting to estimate Qo by Eq. © . As the work func- 
tion of the alloy is unknown, using accordingly 4 and 3.8 
eV for Pb and In, we obtain considerably distinguished 
quantity Qo ~ +13.6 e. However, assuming the exact 
quantity Qo ~ +0.5 e, it is possible to solve the inverse 
problem and to find a work function of the alloy: 3.8012 
eV (instead of 4 eV for Pb). 



B. Structure under voltage V 

Between the emitter (V — 0) and the collector the 
positive voltage V is applied. We consider a central 
electrode-granule in an external electric field. In a weak 
electric field approach we assume, that the ionic sub- 
system of a granule is not deformed, and the electronic 
"cloud" , generated by the own valence electrons, is de- 
formed only. 

The total energy of a granule is the functional of non- 
homogeneous electron concentration, _E[n(r)]. The func- 
tional contains a contribution responsible for the inter- 
action of electrons and ions with an external field, 

e J [n(r) -n,(r)](E -r)d 3 r. (7) 

For simplicity, we suppose that the ion distribution rij (r) 
is spherically symmetric. 

Let's write down an electron distribution of a granule 

as 

n(r) = n Q (r) + 5n\{r) + 5n 2 (r). (8) 

Here, no(r) is the electron density of neutral cluster in 
the absence of the external field, 

J no(r) d 3 r = N , 

Sni is the perturbation arising from the charging of gran- 
ule, and 5ri2 (r) is the next perturbation arising from the 
external field which responses for the polarization of a 
neutral granule, 

J 5 ni {r) d 3 r = AN, J Sn 2 {r) d 3 r = 0, (9) 

where AA^ > and AJV < correspond to negatively and 
positively charged granule, respectively. We assume, that 
functions no(r) and ni(r) are spherically symmetrical, 
and 712(1") is axially symmetrical. Then one can expand 
the -E[n(r)] in the functional Taylor series down to the 
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second order of smallness with respect to 8ni and 6ri2, 
E[n(r)] = E[n (r)\ + ]T J 5n d 3 r+ 

\l a , S lf t ,J n J (r)5n k (v')d\d i r' + ... 
4^ J J Snj(r)5n k {r') 

(10) 

Here the functional derivatives are taken at n(r) = no(r), 
and indexes j and k runs 1 and 2 according to the defini- 
tion (JSJ). The zeroth-order expression E[n (r)] = E 00 is 
a total energy of a cluster before the charging (AN = 0) 
and in the absence of the external field (E = 0). The 
functional derivative 

8E/8n(r) = /x + e(E • r), (11) 

where the constant \x is the chemical potential of electrons 
of a granule. For large clusters, — fi — W in the absence 
of charging and the external field. 

Finally, in the semiclassical approximation, we get 

E = E aQ + fiAN-eANrjV+(AN) 2 Ec/2-aE 2 /2 (12) 

(see Appendix B). 

Solving separately the electrostatic problem for the 
same structure with fraction of the voltage V > (Fig. 
P), we obtain 

e 2 e 3 (d c + e 1 L/2e 2 ) + , . 

Here, L = 2R, H for a sphere of radius R and a disk of 
thickness H, respectively, e\ = e, t% = 63 = 1. Under Eq. 
(|13fl one can find the values r) + < 0.65 and rj + < 0.55 
in experiments d H H 0, @ for the spherical-like and 
disk-shaped clusters, accordingly. 

Now we examine the problem of critical surplus charges 
of a cluster in the presence of an external voltage. For 
convenience, further we use a contraction n = AN. 

C. Coulomb instability of a cluster in electric field 

It is necessary to note, that even the vanishing external 
electric field leads to the instability of a cluster because 
of the possibility of tunnelling of electrons. We assume, 
that a cluster relaxed in a metastable state over a period 
of time which is much smaller, than that between acts of 
tunnelling. As a result of the charging, the intrinsic me- 
chanical "stress" leads to the "Coulomb explosion" . This 
problem was described in Ref. 13] for a single spherical 
cluster in absence of the field. Extending these results, 
one can write the following expression 

T- {(-fJ-e,i + \er)V\)R/e + e/2} (14) 

for the critical electronic or ionic charge in quasi-classical 
approximation. For the range V =(0, 2V) we have: 



i) 77 -C 1. Transitions of electrons between the emit- 
ter and the cluster occur more often, than between the 
cluster and the collector, therefore the electrons are ac- 
cumulated on the cluster. In this case their maximal 
number is 

n max ~ W e R/e 2 + 1/2, 

where W e = W e o — fi e i/R and n max ~ +2.5 — h6.5. 

ii) 77 1. Transitions of electrons between the clus- 
ter and the collector occur more often, than between the 
cluster and the emitter, therefore on the cluster the defi- 
ciency of electrons is observed. Using the ion work func- 
tion, this number determines as 

nmin = -(V^ + |e77U|)i?/e 2 -l/2, 

where Wi = Wig— [in/R and rt m i n « —4 11. Similarly, 

for V =(-2, 0V) we have: 

i) n min ~ -Wi R/e 2 -1/2 = -3.8 - -10.6. 

ii) u max = (W e + \e V V\)R/e 2 + 1/2 w +3 - +8. 
Below, the whole numbers [n max ] and [n m i n ] bound the 

summation in (|27|l . The effect of spectrum quantization 
can change these numbers no more than in ±1 according 
to the first inequality in Q (see Ref. [T^'l. 

Effective collision frequency of excited electrons in a 
cluster is defined as 

T- 1 = t- 1 + v F R-\ v F = (ft/mr s )(97r/4) 1/3 , (15) 

where r is a relaxation time in the bulk of the metal, 
caused by electron-electron collisions (r x 10 14 = 6.23 s 
for Au at T — 75 K 19]), and v F is the electron velocity at 
the Fermi surface in the bulk. The estimation performed 
in Ref. [2(| gives a preferred electron collision on walls 
of a dot, therefore r e ~ R/up. It leads to r e Ae ~ (0.52 — 
0.17) h, i.e. to a broadening of levels. 

During the resonant tunneling, the discreteness of the 
spectrum can be revealed only at low temperature (sec- 
ond of inequalities (|TJ). 

The electron thermalization occurs much faster than 
acts of tunnelling. "New " electrons fill up a number of 
own electrons, changing their distribution and, accord- 
ingly, the chemical potential. This state of the cluster 
will be a starting state for the next act of tunneling. 

III. BASIC ENERGY AND KINETIC 
RELATIONS 

We assume, that the total energy of all three electrodes 
E does not change during the tunneling. In the case of 
transition of 5N electrons from the emitter to the granule 
(containing n "surplus" electrons), we have from i|12|) 

SE = -SN ~? + 6Ns p 
(-e) 2 

+ -hJr K« + ^O 2 - « 2 ] - e6Nr] + V = 0. (16) 
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In this equation we take into account that the electron 
ionizes from the level e e on the emitter (whose capaci- 
tance is equal to zero) and then sticks to the level e p in 
a granule with capacitance C. 

By analogy with Ref. |2lj|, using Eqs. ljlU|l and ijT^t for 
5N = 1 , and than taking into account a contact potential 
difference (0, for emitter-granule transition we have 



e e = e p + E c (n + 1/2) - e V +V, 



(17) 



where e p = e p — e5(j). The arrow on the top indicates the 
energies which are determined by transfers according to 
Fig. 1. We suppose, that n = n(V) and n = at V = 0. 
However, the granule is charged by a Qq before voltage 
applied. Therefore, we assume that n is the result of the 
applied voltage only. 

For granule— emitter transition we have 

P = e p + E c (n -1/2)- erj+V. (18) 

Similarly, for the granule— collector and 
collector— granule transitions we have 

P = e p + E c (nT 1/2) + e{l- V +)V. (19) 

Here the upper/under arrows at the left correspond to 
the following signs on the right. Independently of n the 
relation 

Te *Te 771 c 7c 

£ — £ — £jQ — £ _ £ 

takes place. It agrees with well-known quasiclassical re- 
lation for clusters, IP - EA = E c (see, e.g. Ref. [13}"). 
Eqs. I|17|) _ i|19fl represent a "gold rule" for transitions. 

The tunneling of a single electron through barriers is 
determined by the tunnel rates r e,c which depend on the 
junction geometry and the voltage fraction r\. In general, 
their evaluation is far from a trivial problem 0, 101. We 
assume that they are small and the temperature is not 
too low, i.e. 

k B T > h{T e + T c ) < min{A£ p , E c }. (20) 

By analogy with the theory of Ref. 0] , we introduce 
the partial tunneling rates from electrodes to a granule 



2^r(?)/(? + ^)[i-/(?-^)], (21) 



< = 2^r(?)/(? + ^)[l-/(?-^)], (22) 

p 

and from a granule to the electrodes 

< = 2^r(?)[l-/(? + ^)]/(?-^), (23) 

p 

< = 2^r(£ t )[l-/(£ t + T^{>)]/(?- / i), (24) 

p 



where the factor 2 takes into account the spin degenera- 
tion of levels in electrodes. In view of the applied voltage 
(and charging for a granule) the spectrums (see Eqs. l(T7|) 
- (Iiyfl ) and the chemical potentials are shifted in distri- 
butions © and 



He = M 



E c (nTl/2)-e V + V, 



fi c c = ti-e84>+E c (n±l/2)+e(l-r) + )V, Wy = W^+eV. 

As the first approximation of the perturbation theory 
|13| . for small V, fJ, is determined not only by the formal 
shift of the well depth, but also by the number of conduc- 
tion electrons in the state (N = No+n q , n q = n+[Qo]/e). 
The use of the chemical potentials is correct in a quasi- 
equilibrium state, i.e. when the intervals between acts 
of tunneling are much longer than the relaxation time 
of a granule. It is also supposed, that the external field 
and the Coulomb blockade do not remove degeneration 
of levels. In this case, at T — 0, probabilities l(2*T|) - l|2^|l 
will be nonzero for the intervals: 



& < £ e < —Wy 



< 



< 



< 



< -W$ <e c < n c c . 



Let's denote the total electron transition rates to a 
granule and back on electrodes, as 

<* = < + <> < ut =<+<- 

In the limit of weak tunneling, the probability P n of the 
finding of n for the above mentioned electrons at central 
electrode is defined by the master equation ^(| 

P n = <£ t P n+1 + <_! P n _! - (< + < ut ) P n . (25) 

The requirement of the stationarity, P n = 0, gives the 
recurrent relation 



P„ +1 =P n </<£ 



(26) 



The dc current flowing through a quantum granule 
(with restriction on its instability IjMf) ). is determined 

as 



"max>0 



w, 



nmin<0 



"max>0 

-e ^2 P> I 

"min<0 



(27) 

Let's consider an exotic case of the "strong quantiza- 
tion" for electron spectrum [l6}: 

Ae p > E c - 

This regime is hypothetically reached by a significant in- 
crease of the cluster capacitance (the cluster shape must 
be changed to the needle-like or disk-like one under the 



condition that its volume is fixed (see, e.g. Ref. [T3, ex- 
periments with disks |7fl and nanowires [22j])). Thus the 
residual charge Qo in which provides a contact poten- 
tial difference, is proportional to the capacitance and can 
have large magnitude. When the voltage is applied, the 
charge, which is caused by the transferring surplus elec- 
trons, is much less than Qo- Therefore it insignificantly 
influences the cluster energetics. In reality, the inequality 
Ae p » Ec is not possible for the atomic chain [y| ■ Nev- 
ertheless, this case is useful from the methodical point of 
view to analyze the current gap of J — V characteristics. 

As an assumption, we use the fixed tunnel rates at the 
Fermi level in the emitter. It is correct for the small 
voltages, eV <C W. Neglecting in JUJ - (OJJ terms ~ 
Ec, it is easy to obtain the result, similar to Ref. 16]: 

I = I Q J2[.f(e e + W^-f(e c + W^], (28) 




where I Q = 2er e r7(r e +r c ). At T — ► 0, it is convenient 
to write Eq. I(28jl as a "combination" of the staircase 
functions: 



I/I Q — » ]T m p + Wo) - 6(( p + W$)] , 



(29) 



where for shifted spectrums in the emitter and the collec- 
tor the folowing notations are introduced: £ p = e p — er) + V 

and Cp = e p + e (l - V + ) v - 

Let's remind, that the expressions in this section are 
written down for V > 0. In the case V < 0, the I(V) can 
be received, if we set V — on a collector and V > on 
the emitter. Now for Eas.l(T3|). JTTJ-lUnj) we have 



n 



e = e. 



£2£3(eijc/e3 + eiL/2e 2 ) 
eie 2 d c + £ie 3 L + e 2 e 3 d e '' 



E c (n± 1/2) -e(l-r]-)V, 



(30) 



e e =e p + E c (n =p 1/2) + eufV. 

Also using 

W£ = IU C , Jfc = ti- eScp + E c (n ± 1/2) - e(l - rf)V, 

^ = n - e5<j} + E c (n T 1/2) + erfV, W{, = + eV 

and I — > —J, for derived dependence I(V) at V > it is 
necessary to make mirror reflection relative to V — on 
area V < 0. In such a case, for example, value of the Iq 
in Eq. (|29|l on V > differs from Iq ox\ V < Q since the 
tunnel rates are different. 

In the general case, for calculation of I — V 127(1 it is 
necessary to know probabilities P n . Their statistical de- 
termination is a complicated problem |23| . In the exper- 
iments, the size of the cluster and its location are known 
only approximately, therefore detailed calculations of P n 
are not suitable. Using the recurrent relations we can 
find the ratios Pwo/Pq. 



FIG. 3: The calculated size dependence of the residual charge 
Qo © for the structure based on the clusters of various shape: 
sphere (dotted line) and disk (solid line). For illustration, Qo 
of magic sphere AU439 is marked as x . 



IV. APPLICATION AND DISCUSSION 

The cluster is charged positively before the application 
of voltage. The size dependence of a charge Qo(N ) for 
referred gold clusters is demonstrated on Fig. For 
the above mentioned sizes of spherical clusters, Qo < e. 
However, Qo can accept a values larger then e for the 
disks of monatomic thickness. Additional charging of the 
cluster can lead to the Coulomb instability, because the 
quantity Qo is close to a critical charge 

Moreover, cluster's anomalous electrostriction is pos- 
sible as a result of the charging 24] . 

Setting the collector— granule distance d c , parameter 
(3 = r e /r c and using the recurrent relation (|26|l for Eq. 
(|27p. it is possible to calculate the reduced dc current / = 
I/(ePoT e ). We do not evaluate separately the threshold 
voltages, in our scheme it appears automatically. 

The results of calculations of the I — V characteris- 
tics for the structures Au/Au_/v /Au, based on spherical 
clusters, are presented in Fig. For completeness of 
analysis, the voltage behavior of the reduced probabili- 
ties P n (V) = P n /Po and the stream Aw„ = w e n — are 
given also. 

The current jumps are stipulated by the jumps of 
P_i(V) and Au>o(V), because the current is formed by 
their product. Making use of the equality I = J2 n In(V) 
in accordance with Eq. I|27|l one can fix also the "thresh- 
old" values n. 

As one can surmise, the jump of probability P_i(V) 
causes the current jump in the threshold voltage Vq+. 
Note that the jumps of the stream Awo at V = Vb+ 
and Vo- are determined by jumps of P_i(V) or Awo, 
respectively. 
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FIG. 4: The current -voltage curves (solid lines) and its com- 
ponents, calculated from Eq. I^TJ (/3 = 1, rj + = 0.1, T = 30 
K). Awn(V) is given in F e units. 



As is seen from Fig. 4, the role of partial current com- 
ponents (with \n\ > 1) grow with increasing Nq. The 
charging leads to energy shift of spectrum according to 
Eqs. (|17f) - H19|) . Thus the different parts of a spectrum 
are involved during tunneling. The electron chemical po- 
tential of magic cluster AU40 does not coincide with a 
energy level at a zero voltage. 

The current gap width AV g — Vb+ + \Vo-\ in all cases 
is determined by values n = 0,-1, and the boundaries 
can be defined by the position of lowest unoccupied level 
as 



l^o±| 



Ec/ 2 + Ae 



(2-r?±) ' 

where Ae = {{j, p — e LU and 0} for the magic and non- 
magic clusters, accordingly. The probability P_i prevails 
over P+i, because the "granule— collector" stream is more 
than a "emitter— granule" one and the granule is charged 
positively (i.e. n < 0). 

Within the applied voltage the I — V characteristics 
versus -q + are shifted to the right and the gap width 
decreases a little. The calculated I — V curves of the 
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FIG. 5: Calculated I — V curves at T = 30 K for structure 
based on spherical clusters. For presentation the curves are 
shifted slightly on a vertical. 



structure Au/Augoo/Au for fixed rj + and different j3 are 
shown in Fig. I n this model, the current gap is prac- 
tically independent on /3, however, the current jumps are 
strongly dependent on the value of /?, which, in its turn, 
has no influence on threshold voltages. 

In order to illustrate our results, in Fig. we com- 
pare the size dependences AV^(Ao) calculated from Eqs. 
(|27Jl and Ij28(l for spheres and disks. The largest quan- 
tities AV g correspond to the magic granules, for which 
Ae 7^ 0. For the case of "strong quantization" the size 
of a current gap for non-magic clusters is equal to zero 
explicitly, because the emitter Fermi level is in line to 
the closed levels in cluster. Calculations demonstrate the 
non-monotonic dependence AV^(7y). These results shows 
also, that a charging leads to the growth of a gap. 

The actual forms of current gap for the structure based 
on magic disk Aui7g (R ~ 3.5 nm) are plotted in Fig. [7| 
Our calculations show that the dependence AV^nV) is 
slightly non-monotonic. However, in experiments 0, Q 
the gap varied as 0.8 — » 0.4 — > 0.7 V for the cyclic vari- 
ation ofd c «1^2^lA. The reasons of such con- 
siderable difference, apparently consist in the following 
effects: broadenings of levels, amplifications nonlineari- 
ties in the strong electric field, and in energy dependence 
of tunneling rates. At high rates the capacitance ceases 
to be classical and can strongly grow (Ec — > 0) 0, |2(| , 
showing non-monotonic dependence from T c . It means, 
that in a reality we deal with the intermediate cases (be- 
tween limiting estimations from Eqs. (|27|) and l|28l) in 
Fig. 7). Only in these situations the observed behavior 
of gap width is more or less clear. 

Let's discuss other features of the tunnel construction. 
In spite of the fact that the emitter and a collector are 
made of one material, the chemical potentials of electrons 




r,(y) 

FIG. 6: The current gap calculated from Eq. 12711 (d c — 
2 Aand f3 = 10). Solid lines show the gaps calculated from 
Eq. I|28^ for the case of "strong quantization" . 



are not equal to each other: the emitter is represented 
by a thick film of Au (111), and a collector is a poly crys- 
tal of Au. Their work functions are different [27]]. We 
assume also that W < Wo for low-dimensional systems 
(l3l | . Except for it the emitter is covered with a dielectric 
film, that also influences a electron work function. We 
can estimate this contribution. 

Proceeding from indirect experimental measurements 
|28j . the work function decreases with growth dielectric 
constant e of coating. The calculations of the electron 
work function Wd for cylindrical nanowires in a dielec- 
tric confinement are done in Ref. 29): W4 decreases 
approximately on 20% at magnitude as e rise from 1 to 
4. The basic contribution thus can be related to the 
change of electrostatic dipole barrier which contribution 
to a work function of system gold— vacuum makes up to 
30% |2^ . Hence, this contribution also makes a upper 
limit of the change of Wd for metal— dielectric— vacuum 
system. Owing to Wd < Wq the inequality W > Wd is 
possible, that leads to negative charging of the cluster 
before the application of voltage. The film— cluster con- 
tact also will change the cluster energetics. At last, the 
metal-nonmetal transition for gold cluster can appears 
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V. SUMMARY 

In this paper, we presented an approach for the calcu- 
lation for the I — V characteristics of SET based on the 
metal clusters. 

In the framework of the particle- in-a-box model for the 
spherical and disk-shaped gold clusters, the electron spec- 
trum was calculated. In this model, the work function of 
clusters is smaller that of semi-infinite metal electrodes. 
It resulted in the appearance of a contact potential differ- 
ence between a cluster and electrodes. Residual charge 
is equal to non- integer elementary charge e, which corre- 
sponds to the charge of cluster in chemisorption regime. 
For the small spherical clusters, positive charge is less 
than e. However, this charge can accept a values larger 
than e for the disks of monoatomic thickness. Additional 
charging of the clusters can lead to the Coulomb instabil- 
ity, because it is close to a critical charge. The charging 
results in the energy shift of the spectrum. 

The current-voltage characteristics were analyzed tak- 
ing into account the contact potential difference. The ap- 
proach was applied to calculate the trapped offset charges 
that determine the transport behavior of molecular-like 
structures. For single-electron transistors based on the 
small gold clusters the current gap and its voltage asym- 
metry were computed. The largest quantities gap cor- 
respond to the magic granules. For the case of "strong 
quantization" , the size of a current gap for non-magic 
clusters is equal to zero explicitly, because the emitter 
Fermi level is in line to the closed levels in cluster. The re- 
sults shows, that a charging leads to the growth of a gap. 
The evaluations demonstrated a slightly non-monotonic 
dependence (in comparison with the experiments) of the 
gap size from cluster - collector distance. At high tunnel- 
ing rates the capacitance ceases to be classical and must 
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strongly change. 
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In addition to the spin degeneration, there is a double 
degeneration with respect to the sign of index m, since 
k n ,m = k n ,-m- Further, the spectrum of cluster is de- 
noted 1, 2, 3... is the number of an one- electron 
state. All levels are numbered in order to increase ener- 
gies. 



APPENDIX B: ENERGY OF CLUSTER IN 
EXTERNAL ELECTRICAL FIELD 



APPENDIX A: ELECTRON SPECTRUM IN 
CYLINDRICAL-LIKE CLUSTERS 

As an approximation, the profile of the one-electron 
effective potential in the cluster can be represented as a 
potential well the depth Uq < 0. The three-dimensional 
Schrodinger equation for a quantum box can be separated 
to the one-dimensional equations. The spectrum of wave 
numbers in a spherical and cylindrical potential wells are 
determined from the continuity condition of a logarithmic 
derivative of the wave function on the boundaries. 

For a disk of radius R and thickness H it is necessary 
to solve the equation: 



Irn (k nm R) 



K' m (>c nm R) 

■^m ( l^ri m R ) 



(Al) 



Here I m is the Bessel function, K m is the McDonald func- 
tion, the stroke denotes a derivative over an argument, 
k nm = yKQ — Xnrm foko = \J 2m e \Uo\, and m e is the elec- 
tron mass. The n = 1, 2, 3... number the roots of the Eq. 
fAlj) for the fixed m = 0, ±1, ±2.... 

Using the perturbation theory , the finite well prob- 
lem can be reduced to the infinite well one: 

knm — fc^m + k nn \ + \k nm /k nn \\ <C 1, 

where k n °m defines the spectrum of infinitely depth well. 
Numbers k n ° m are determined by solutions of the equation 



In 



R) = 0. 



In the first approximation we have 



R^lK> n {x n ZR) 



(<>) 



The estimation gives kn m < 0.07fcim, confirming suffi- 
cient accuracy of the theory. 

Quantization of the wave vector k s along the cylinder 
axis is determined by the solution of the equation: 

k s H — sir — 2 arcsin(fc s //co), 

where s is the integer number. Neglecting the area near 
cylinder edges, the energy spectrum is calculated by a 
simple way as follows 



Uo + ^(k nm + k 2 s ). 



Using spherical coordinates, we remove the center 
point z = in a center of a granule, and we direct a 
z axis from a collector to the emitter under the conser- 
vation the potential difference between them. Then an 
electric field E = Ez, where z is a unit vector along an z 
axis. 

As the surplus charge is effectively distributed over a 
surface, it is quite reasonable to use 



5m(r) = A5n(r - R), 



(Bl) 



where 5b {f — R) is the Dirac delta— function. In spite 
of the fact that the form i|BID corresponds to the total 
screening of a surplus charge inside a granule, this form 
is rather convenient for the calculations. Then, we use 
the linear response approach (see, e.g. Ref. ^j) 



5n 2 (r,6) = Y(r)Ecos6. 



(B2) 



The constant A in Eq. (|B1|I and spherically symmetric 
function Y{r) in (|B2p are determined from the normal- 
ization condition and a global minimum of the func- 
tional, 5E[n(r)] -> 0. 

One of the terms, interesting for us, is 

— e J 5ni(r)ip(z) d 3 r, 

where <p is an external electrostatic potential. In the case 
of V > and vacuum collector-emitter space <p(z) = 
V(z — d e — L/2)/d, where d = d e + L + d c . 

After the integration in spherical coordinates, the 
term, which is proportional to z, vanishes, and, as a 
result, we have — eAN-qV, rj is a fraction of a voltage. 
Other three terms 



<5ni(r)(5n 2 (r') + 5ni(r)5ni(r') + 5n 2 {r)5n 2 {r') 



d A rd 



■3„/ 



(B3) 

give a basic contribution to the second order of expansion 
(Unf . The first integral in Eq. I|B3}1 for the functions I|B1|) 
and l|B2|) vanishes after the integration on corners, second 
and the third ones were calculated earlier at the definition 
of the ionization potential (see, e.g. Ref. and polar- 
izabilities of a cluster a = — (4tt/3) J^°Y(r)r 3 dr = i?^ ff 
[l7j . Finally we obtain Eq. (tT2l . 
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